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This research examines an admissible search algorithm-dependent A-star 
strategy and takes into account redundancy and cost/benefit analysis under 
normal operating conditions. The goal is to allocate a phasor measurement 
unit (PMU) for maximal observability of the interconnected power network. 
To determine the fewest number of PMU required to make the connected 
power network totally observable using redundancy analysis, the A-star 
approach is utilized. The redundancy analysis of the power network is 
carried out in order to determine the appropriate PMU placement, which 
results in the acquisition of total power network observability and reliability. 
To put the suggested technique through its paces, it has been tested on IEEE- 
standard test systems such as IEEE-14 bus, IEEE-30 bus, New England-39 
bus, IEEE-57 bus, and JEEE-118 bus. The results obtained using the 
suggested methodology are compared to those obtained through standard 
literature research. The experimental findings of the suggested method 
revealed the resilience and accuracy of the A-star algorithm as well as its 
effectiveness in achieving maximum observability of the connected power 
network. 
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NOMENCLEATURE 

PMU : Phasor measurement unit ILP Integer linear programming 
TEEE : Institute of Electrical and Electronics Engineers ZIBs  : Zero injection buses 
SCADA : Supervisory control and data acquisitions RBs : Radial buses 

OPP : Optimal PMU placement RA : Redundancy analysis 
RCOs : Random component outages CBA : Cost/benefit analysis 
MCDM : Multi-criteria decision-making MR : Maximum redundancy 
BPSO : Binary particle swarm optimization 


1. INTRODUCTION 


Phasor measurement unit (PMU) is an online monitoring instrument in the wide-area measurement 
system, which is used for monitoring and supervisory purposes in leading and complex interconnected power 
networks. Till now, supervisory control and data acquisition (SCADA) is used for state estimation, but it 
provides asynchronized measurement dominant towards the unreliable analysis of power network 
conditions [1], [2]. In addition, the scan rate of data is 2-4 samples per cycle which makes the SCADA system 
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incapable of measuring dynamic/transient disturbances in interconnected power networks. PMU first came into 
existence in the mid-80s, which utilizes synchronized signals and provides synchronized measurements for 
actual time phasors of voltage and current. Synchronization is attained by the actual time sampling of voltage 
and current waveforms using timing signals from the global positioning system (GPS) clock [3], [4]. The 
observability of the interconnected power network is attained by the real-time data which is taken from PMU. 
PMU located at each bus in a connected power network efficiently observes all the states. But it is unwise and 
uneconomic to place the PMU at each bus in a power network as PMU and its communication services are 
expensive. Hence, an appropriate methodology is required for the optimal PMU placement problem (OPP) in 
the interconnected power network. The applicability of the present paper is restricted to the use of PMU for the 
complete observability of connected power networks. A connected power network is completely observable 
only when all of its conditions are evaluated separately, as given in [5], [6]. 

The appropriate site selection for placement of PMU in interconnected power networks has become a 
vigorous research issue [7], [8]. In recent years, various classical soft computing techniques and approaches 
have been published in the open literature for the optimal placement of PMU [9]. The authors propose a novel 
binary search algorithm-based methodology in [10] for the OPP problem for the complete observability of the 
interconnected power network states. Also, redundancy analysis is evaluated for the appropriate location of 
PMU. A novel three-step methodology for site selection of PMU using network connectivity data is proposed 
in [11]. The procedure considers PMU at each bus in a connected power network for observability of the states. 
Step 1 and Step 2 of the methodology decide: i) the removal of inappropriate bus locations of PMUs and 
ii) essential bus locations where PMUs are kept. Step 3 minimizes the number of PMUs using a pruning 
application. Abiri et al. [12], the ramification of the channel capacity of PMU on their site selection for total 
state estimation of the power network is suggested. Initially, the standard assessment of the grid network is 
formulated and then the revised methodology is indicated for the OPP problem. A genetic algorithm (GA) based 
approach is proposed in [13] for the OPP problem. It determines the minimum number and optimal locations of 
PMU and their geographical distribution for the full observability of the power network. The author proposed an 
intelligent search-based approach for the OPP problem in [14]. Also, redundancy measurement is taken into 
consideration for minimizing the number of PMU. The authors proposed a novel methodology for the site 
selection of PMU suffering from random component outages (RCOs) in interconnected power networks [15]. 
Using the RCO model, the optimal PMU locations are preferred to minimize the state estimation error 
covariance. 

Gómez and Ríos [16], the authors suggested two conditions for OPP issues. Firstly, an integer linear 
programming (ILP) based model for the OPP problem is presented. This methodology determines the minimum 
number and locations of PMU at each step while exaggerating the interconnected power network observability 
at each time. Secondly, a procedure is proposed to find the number of significant buses to be analyzed for 
dynamic/transient stability monitoring. An original method based on integer linear programming (ILP) and 
multi-criteria decision-making (MCDM) is provided in [17] to ensure the observability of power networks even 
during branch outages or PMU collapse. For the optimal location of PMU, ILP is used and then PMU location 
is prioritized by the proposed MCDM model. The authors propose a new investment decision model in [18] for 
the placement of PMU to ensure a complete estimation of the interconnected power network. Ramachandran 
and Bellarmine [19], a new meta-heuristic-based approach, namely the fruit fly optimization technique, is used 
to find the minimum number and location of PMU for the entire observability of the interconnected power 
network. A simple methodology for OPP problems in connected power networks using ILP is suggested [20]. 
A combination of graph theory and GA is used for the OPP problem to ensure state estimation of the power 
network [21]. Also, Theodorakatos [22], numerical observability formulation is used to determine the strategic 
locations of PMUs considering branch-and-bound and a binary-coded genetic algorithm. Further, a two-phase 
branch and bound are proposed in [23] to unravel the OPP problem. Later, a modified branch and bound are 
developed in [24] in order to solve the PMU placement problem. 

Kavasseri and Srinivasan [25] deal with the joint OPP problem and power flow measurements to 
provide observability under faulted states in interconnected power networks. Initially, a non-linear integer 
programming technique is adopted to solve the problem, and then it is changed into an equivalent ILP problem 
through Boolean implications. Similar work is done in [26], in which both the stochastic and deterministic 
methods are used to find the fault location observability problem. Liao et al. [27], a hybrid two-phase approach 
is proposed for the OPP problem for a grid-connected power network. In the first phase, candidate locations of 
PMUs are found by using a graph theory approach. Then the local search heuristic approach is developed for 
finding the minimum number and optimal locations of PMU in the connected power network. A two-stage OPP 
method is proposed in [28]. In stage-1, the minimum number and locations of PMU to make grid network fully 
topologically observable and in stage-2, it is tested whether the obtained results of PMU locations lead to 
complete ranked measurement jacobian or not. 
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Xia et al. [29], an alternative approach for OPP with redundancy analysis for complete state estimation 
of the grid network conditions is proposed. A multi-objective biogeography-based optimization technique for 
OPP issues ensures that the interconnected power grid is fully observable [30]. In [31] and [32], the authors 
offered an ILP algorithm for an OPP problem to ensure complete observability of the connected power network 
states using tomlab optimization. Yazdel and Esmaili [33], a reliability-based approach is suggested for OPP 
problems and flow measurements. A GA-based method is used for the OPP issue in [34] and bus ranking 
formulation is developed for finding the minimal number of PMU necessary to make interconnected power 
networks entirely observable. A BPSO is recommended to solve the OPP problem [35], [36]. 

The purpose behind the whole research work is to propose a new method for solving the OPP problem. 
The researchers have adopted several techniques to solve the OPP problem, as stated in the literature survey. 
However, most of them are based on network topology and the evolutionary approach. The present paper has 
utilized the concept of admissible searching algorithm-based A-star approach for the OPP problem in order to 
determine the locations of PMUs, considering redundancy and optimally cost analysis. The main contributions 
of this present task are: 

— To perceive the minimum number and optimal locations of PMUs to make the interconnected power 
network topologically observable. 

— OPP methodology stated in this paper assures that the power system is entirely observable under normal 
operating states. In this work, the authors propose an A-star-based approach to find optimal locations of 
PMU for observability of power network states during normal operating conditions. 

— Redundancy and cost analysis of PMU locations is attained by using the A-star technique in which the 
preferred final solution contains optimal locations of PMU for complete observability of the network states. 

— The proposed approach is tested on IEEE 14-bus [37], IEEE 30-bus [37], New England 39-bus [38], IEEE 
57-bus [37], IEEE 118-bus [37] test systems, and obtained results are compared with the standard referred 
journals published. 

The structure of the paper is as follows: PMU placement rule is described in section 2. Section 3 
explains the implementation of proposed approach to unravel OPP problem for the complete observability of the 
interconnected power network. In section 4, a discussion on case studies and results is performed. Section 5 
concludes the paper. 


2. PMU PLACEMENT RULES FOR INTERCONNECTED POWER NETWORK 
The proposed OPP approach is applied to standard test systems so that PMUs can be placed at 
appropriate locations to estimate the connected power network states. 


2.1. Methodology for optimal PMU placement (OPP) problem 
The OPP issue can be prepared methodically as: 


min Yat CeXk (1) 
s.t.a(X) =b (2) 


Where, N is the number of buses in connected power network for placement of PMUs, cx is the weight factor 
of estimation to the cost of placed PMU at k*” bus, and X is the binary decision variable vector having 
element x, which decides attainability of PMU on k“ bus. Binary decision variable vector is explained by 
in (3) and a(X) is the observability constraint specific whose appearances are non-zero if the analogous 
different bus voltages are noticeable with respect to the given sets of measurements according to the above- 
mentioned rule or its appearances are zero otherwise. 
_ (1if PMU is required at kt” bus 
ga | 6) 
0 otherwise 


The entries in a are defined as: 


Lifk =i (4) 


1, if bus kis connected to bus i 
ani = 
0, otherwise 


And b is a unit vector, set as: 
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2.2. Observability analysis based on PMU placement 

When PMU is installed at a bus, it measures the voltage phasor at that bus and the current phasors of 
all the branches linked to that bus relying upon the number of channels. It is presumed in this work that PMU 
with an adequate number of channels is set up at a bus so that the current phasors of all the branches linked to 
that bus can be evaluated and voltage phasor of that bus too. The voltage phasors at the buses acquainting to 
the PMU setup bus can be resolved using the measured branch current phasors, bus voltage phasor and 
known line parameters [10]. Thus, with the optimal number of PMUs at a subset of power network buses, the 
complete estimation of the interconnected power network is attained. Figure 1 shows the observability of the 
connected power network with the help of PMU. As shown in Figure 1(a), if PMU is placed at B3, then it 
observes not only that bus but also all the acquainting buses, i.e., Bı, B2 and B3, are connected to it. Assume if 
buses Bz, B4 and including ZIB, B3, are observable by PMU, but bus Bı is unobservable, then Bı is observed 
by using (1) as shown in Figure 1(b). Suppose, if buses Bı, B2 and B4 are observed by PMU but B3, a zero- 
injection bus is unobservable then B3 is measured using (1) as depicted in Figure 1(c). As shown in 
Figure 1(d), assume if B3 is a zero-injection bus, then PMU must be placed at ZIB to observe the B3 and its 
acquaintance buses B1, Bz and B4. 


OBSERVED OBSERVED 
BY PMU BY PMU 


PMU 


IMPERCEPTIBLE ^, 
BUSB1IS 
OBSERVED BY 
KCL 


BUSES 
INCLUDING ZIB 
ARE OBSERVED,’ 


` BY PMU |” 
OBSERVED See ai 
BY PMU 
(a) (b) 
OBSERVED OBSERVED OBSERVED OBSERVED 
BY PMU BY PMU BY PMU BY PMU 
H ` G) 
‘ t; 
: ; ZIB 
IMPERCEPTIBLE \, / CONNECTED 
BUSB3IS “| d WITH PMU 
OBSERVED BY 
KCL OBSERVED 
OBSERVED © BI mu 


BY PMU 


(c) (d) 


Figure 1. Observability of interconnected power networks with the help of PMU: (a) PMU is placed at Bus 3, 
(b) ZIB observed by PMU, (c) ZIB is observed via KCL, and (d) PMU is placed at ZIB 


In Figure 2(a), if PMU is placed at bus Bo, then it observes Bı, B2 (connected with PMU), B3 and B4 
but Bs, a radial bus (RB), is unobservable. Then, Bs is measured using (1). As depicted in Figure 2(b), if PMU 
is placed at a radial bus, then it observes only one bus B4 and itself. Therefore, the radial bus is not included 
in the placement of PMU. 
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Po aeeoe 7” TWO BUSES ARE 
OBSERVED 
AFTER PLACING 
PMU ATRADIAL 
BUS 


(b) 


Figure 2. Radial buses (RB) in the connected power network: (a) PMU is placed at Bus 3 and 
(b) PMU is placed at RB 


3. PMUs FOR COMPLETE OBSERVABILITY OF THE INTERCONNECTED POWER 
NETWORK 

PMU has now become the measurement technique of choice for connected power networks. 
The advent of PMU makes the power network fully observable during transient/dynamic states. 
The interconnected power network consists subset of buses i.e. substations, power generating units, or a 
huddle of loads and transmission lines connecting all the buses. PMU placed on power network buses gives 
the same-time evaluation of interconnected power network variables. This has been expedited with the help 
of GPS technology. A PMU established at a specific bus determines the voltage magnitude and phase angle 
of that bus and the current phasors of all branches connected to that bus rely upon the number of channels in 
an interconnected power network [39]. Accordingly, a power network is fully estimable only when all of its 
buses are observable over direct or indirect measurement. In the present work, the topological observability 
approach is used. Topological observability depends on graph theory methodology, where disoriented graphs 
illustrate the interconnected power network. An interconnected power network can be treated as topologically 
observable if not less than one mensuration tree is of full rank [40]. 


3.1. Modelling of zero injection buses 

In an interconnected power network, some buses are neither connected with load nor with a 
generator, and such buses are termed as zero injection buses (ZIBs). ZIBs are modeled as in [11]. At ZIBs, a 
current is not injected into the connected power network. When ZIBs measurement is not taken into account 
for network observability, then pseudo-measurement is used for the connected power network observability. 
As depicted in Figure 3 buses are named as Bı, B2, B3, and B4 where B3 is represented as a zero-injection bus. 
By applying KCL at B3, it gives the: 


l3 + h3 + l43 = 0 (6) 


Figure 3. Modeling of zero injection bus 


3.2. Proposed approach 

The minimum number and optimal locations of PMUs essential for full observability of the 
interconnected power network is found by using an admissible searching algorithm based on the A-star 
approach [41]. An A-star approach first came into reality by P. E. Hart, N. J. Nilsson, and B. Raphael of 
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Stanford Research Institute in 1968 [42]. It is a fortification of heuristic approaches like best-first search and 
the Dijkstra algorithm [43]. It is used in finding the shortest path. It is a probing-based approach like a greedy 
best-first search where a graph traversing algorithm observes the tiniest cost path from a stated initial node to 
the most promising node (out of one or more desirable goals) [43]. To solve the order in which the search 
visits nodes within the tree, it uses a distance-plus-cost node x (usually denoted by f(x)). 

The distance-plus-cost probing is an addition of two functions: 
g(x) = The actual cost to a most promising node x. 
h(x) = Imprecise cost from node x to goal node x. It is a probing function. The testing function never 
considers overestimated costs, i.e., the actual cost of reaching a target node from node x should be > or = h 
(x). It is called the permissible probing approach. The absolute cost of each node is summed up as: 


f(x) =g) +h) (7) 


A-star only bolsters the node if it appears most promising. It only targets to reach the goal node from the 
initial node and is incapable of getting each node. It is optimal if the probing function is permissible. 


3.2.1. Implementation of the proposed approach 

The key goal of the OPP issues is to render the device utterly measurable with a minimum number 
of PMUs. Thus, the node having higher connectivity of the branches should be favored for installing PMUs. 
This approach is used in the depth-first search (DFS) [44], but results give redundant PMUs also, which is 
unnecessary. The DFS is a graph theoretic procedure-based approach that searches each node in an 
interconnected power network leading to an inadmissible path to arrive at the goal node. The A-star approach 
overcomes this problem. To determine the fewest possible nodes for OPP, the A-star algorithm uses it to 
evaluate each node to decide which one should be augmented next [45]. The fewest possible nodes are 
recognized as those having minimum or maximum scores of an evaluation function. Although searching 
about the slightest optimal cost path, the A-star algorithm can alter its search path from the ongoing search 
path to the fewest possible way. This feature makes the A-star algorithm exceptional over other graph- 
theoretical procedures for OPP issues. Figure 4 shows the flowchart of the proposed A-star algorithm for the 
solution of the OPP problem. 


— 


Assemble binary connectivity 
matrix a(NbxNb) using eqn (5) 


Sy 


Initialize first set of estimation for 
g(x), h(x) and iteration number 
(c_iter) 


— 


Create initially binary connectivity 
matrix for g_x (nxn) and h_x 
(nxn) randomly 


Place PMU randomly either on 
bus B or one of its acquaintance 
bus 


Is goal attained? 


Calculate sum of 
PMUs in each 
row for h(x) 


Choose the best locations 
for PMU (Nemu) in power 
network 


Measure the connected power 
network coverage with installed 
PMUs 


Update g_x and h_x 


Is power network 
fully observable? 


Discard PMU 
randomly 


Create a matrix containing the 
indexes of all the connections 
to the node with PMU 


| 


Sorting the number of buses 
on which PMU are connected 


Place PMU at the bus (B = 1) 
which is connected with max no of 
branches in unobservable region 


Iter = Iter + 1 


Figure 4. Flowchart of the proposed A-star algorithm for the solution of the OPP problem 
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3.2.2. Step-by-step procedure 
The implementation of the A-star approach to obtain the optimal location of PMU for a system 
consisting of ‘N,’ number of buses is given in the following steps: 
— Read bus and line data of standard test system. 
— With the help of given data, obtain binary connectivity matrix a(N, X N,) as in (4). For example, in the 
IEEE 14-bus system, the connectivity matrix is obtained as: 


Bus, 
Busy 
Bus3 
Bus, 
Buss 
Bus, 
BUS7 15 
BuSg pp 
Buso 
Busio 
Busy, 
Busy 
Bus43 
Busy4 


SOC COOCOO COR CORRE 
SOC CO OOOO ORKBRHEEBER 
SOC OCOOO COC ORKLK:,D 
SOC CO OR OH ORBRKBLHO 
SOC COC OOO ORKBRHROFRR, 
SCrRRKFEODOCOHHROGOCO 
SCOCCOORKRHOOHKROGCO 
SOC COOORHKOCCOCCCO 
PODORKROHGCOKRGGCO 
SCOOFRKRHOCOCCOCCSO 
SCOOKRKHKOCOCOHCOCGCCSO 
SCrRRODOCOHCOCCO 
PRRODOCOKRCACGCCSO 
PROSDCOFRCOCCCOCCSO 


— Create an initial binary connectivity matrix randomly and carry out the first set of estimations for g(x), 
h(x), cov and rocov. 

— Set ‘counter(c = c + 1)’. The counter is used to keep track of which row of each of the matrices we are on. 

— Estimate the maximum ratios and save ratios values, r = “ if r = 0, then delete the row so that it 
cannot be chosen again. 

— Decrease the counter, ‘rocov_counter — 1’. 

— Find the column index for each bus connection, ‘c_update’. 

— Determine the number of columns of the cov_find (y = size (cov_find)), matrix for later use of doing 
a search through the cov_find matrix. 

— At each step of the ‘for loop’, take the column index and store 2 in the cov_update matrix. Place 1 to 
represent that a PMU has been placed at kt” node. The 2’s are variables to show that the bus is covered. 
Check for less than optimal as: 


2, if bus i is fully observed by PMU, 
y° =} 1,PMU has been placed at node í where y“ is the observability vector at c** iteration. 
0, otherwise 


— Get the previous heuristic cost h(x) associated with placing the previous PMU 
‘h(x) (r_update, c_update)’. 

— Compute the row and column index of each PMU placed in the last row of the cov matrix (rowgq, colw). 

— Obtain the number of PMU placed, ‘cov_size’. 

— Add a new row to g(x) whose value equals the number of PMUs placed times the number of buses in the 
system. 

— Create a matrix containing the indices of all the connections to the nodes with PMU, 
‘cov_find_additional’. 

— Determine the number of buses on which PMUs are connected, Nemu = size (as in step 14). 
Npmy = {2, delete row and column of the heu_adj that represents a covered node 0, otherwise 
where Npyy is the number of buses on which PMU are connected. 

— Sum up each column of heu_adj, so that it can be used to find the new row of h(x). 


> 0, checks for nodes that still have the connection 
h(x) }= h(x)_update, add the number of connection that node has to be 
the coverage value of the previous decision 


— Make sure not to place same PMUs in a different order, such that h(x) (r_update, c_update) = 0. 
— For complete observability (COpyy) of the connected power network 
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1, checks for full coverage by checking to see if there are any zero nodes 
(COpmy) = in the last row of hy 
0, otherwise 


— Display the PMU locations (Npyy,) 


3.3. Redundancy analysis 
After determining the optimal locations of PMUs, the explanation for redundancy analysis (RA) can 
be presented as: 


RA = EREM" ON ae (8) 


Where Npmy is the total optimal number of PMUs, a is the binary connectivity matrix obtained from (4), Npyu, 


is the PMU locations at interconnected power network buses, which is obtained from step 19 of the proposed 
approach. In (8) allows the redundancy analysis for all the placement sets of optimal locations PMU. 


3.4. Cost/benefit analysis (CBA) 
In this section, the cost analysis associated with the optimal placement of PMU has been 
shown [46]. The cost analysis is formulated by (9). 


CBA = yi," ch Nfmu, + FC (9) 


Where, Npyy and Npyy, are explained in subsection 3.3, c, (=$60,000) is described in subsection 3.1 and 


fixed charge (FC=$400,000) mainly refer to charge for hardware and software facilities. 


4. RESULTS AND DISCUSSION 

The optimization problem, described in the previous sections, is developed and executed in 
MATLAB software on an Intel Core i3 processor-based system with 2.2 GHz clock speed and supported by 
3 GB of RAM. The proposed algorithm for the OPP problem is applied on standard IEEE 14-bus, IEEE 30- 
bus, New England 39-bus, IEEE 57-bus, and IEEE 118-bus test systems. In the present work, B1, B2, B3,...., 
Bn are defined as bus numbers where n = 1, 2, 3,..., N. 


4.1. Test system 1: IEEE 14-bus system 

The IEEE 14-bus test system is taken into consideration as test system 1. The data of the system is 
shown in Table 1. The proposed algorithm is used to determine the OPP problem. Simulations are executed 
by ignoring and considering zero injection buses. Results are presented in Tables 2 and 3. The proposed 
method gives multiple solution sets for the OPP problem, as shown in Tables 2 and 3 ignoring and 
considering ZIBs. The IEEE 14-bus system has three optimal PMU placement sets with four optimal 
locations in each set ignoring ZIBs, i.e {B2, Bz, Bri, Biz}, {B2, Bo, Bs, Bo} and {B2, Bs, B7, Bo}. For more than 
one solution, redundancy analysis is taken into consideration. The maximum value of redundancy analysis is 
taken as the best PMU placement set. Further, the installation cost is also taken into account for the 
placement of PMUs. As shown in Table 2, if PMU is placed at Bz, Bs, B7, and Bo, then the maximum 
redundancy (MR) is 3. 

Here By is observed three times, Bs, B7, and Bo are observed two times and the rest of the buses 
observed once. In the IEEE 14-bus system, the maximum redundancy value for the first set is 16, for the 
second set is 17 and for the third set is 19 ignoring ZIBs. The set having maximum redundancy value is 
considered for PMU placement according to the proposed algorithm. So {B2, Bs, B7, Bo} is the best placement 
set for PMU ignoring ZIBs and {B2, Bs, Bo} is the best location considering ZIBs. Bold letters give the most 
promising sites for the installation of PMUs. 


Table 1. System data of IEEE 14-bus system 
#oflines #ofZIBs #ofRBs Loc. of ZIBs Loc. of RBs 
20 1 1 B, Ba 
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Table 2. Simulation results for the IEEE 14- bus system ignoring ZIBs 
# of PMUs Loc. of PMUs # of times each bus observed MR RA CBA($) Comp. time ins. 


4 By,,By,By, Bı 11121211111111 2 16 640000 0.0394 
4 Bz, Be, Bg, Bo 11122121111111 2 17 640000 0.0429 
4 Bz, Bg, By, Bo 11132121211111 3 19 640000 0.0351 


Table 3. Simulation results for the IEEE 14- bus system considering ZIBs 
#ofPMUs Loc. of PMUs _# of times each bus observed _MR RA _CBA($) Comp. time in s. 


3 By, By1, B43 11111211111111 2 12 580000 0.056 
3 Bz; Be, Bo 11122111111111 2 15 580000 0.091 
3 Bz, Be, Bo 11122111111111 2 15 580000 0.078 


4.2. Test system 2: IEEE 30-bus system 

IEEE 30-bus system is considered as test system 2. The information of the system is shown in Table 4. 
The proposed algorithm is used to determine the OPP problem in the IEEE 30-bus system. In test system 2, 
IEEE 30 bus system has three sets for PMUs placements. The first set is {B,, Bs, Bg, B10, B12, Big, B23, B27}, 
the second set is {B1, Bs, Bs, B12, B17, Bio, B22, B23, B27, B28} and the third set is 
{B,, Bs, Bg, B10, B12, Big, B24, Bos, B27, B28} ignoring ZIBs. The maximum redundancy value for the first, the 
second, and the third set is 34, 39 and 43 respectively as shown in Table 5. Hence, the third set is seemed to 
be for the placement of PMU. But, the set with a minimum number of PMU should be taken into 
consideration for the placement as the installation cost of PMUs is less. Further, the maximum redundancy is 
3. Here one bus observed three times, three buses are observed three times and the remaining buses are 
observed once. In Table 6, eight locations in each set are found considering ZIBs. In this case, the set with 
maximum redundancy value is considered for placement of PMU as installation cost is uniform for each set. 
The most promising locations to place PMUs optimally in the IEEE 30-bus system are shown in bold. 


Table 4. System data of IEEE 30-bus system 
#oflines #of ZIBs _#of RBs Loc. of ZIBs Loc. of RBs 
41 6 3 Bo, Bo, B22, Bas, B27, B28 B11, By 3, B26 


Table 5. Simulation results for the IEEE 30-bus system ignoring ZIBs 


# of PMUs Loc. of PMUs # of times each bus observed MR RA CBA($) Comp. time in s. 
8 B,, Bs, Bg, Bio, Biz, 12111211111111311111 3 34 880000 0.297 
Big, B23, B27 1111111211 
10 B,, Bs, Bg, B42, B17 Bap 12111211131211221111 3 39 1000000 0.301 
B22, B23, B27, Bag 1112112211 
10 B,, Bs, Bg, Bio, B12 Bag 12111312111111211111 3 43 1000000 0.278 
B324 B25, B27, Bag 1212313311 


Table 6. Simulation results for the IEEE 30- bus system considering ZIBs 


# of PMUs Loc. of PMUs # of times each bus observed MR RA CBA($) Comp. time in s. 
8 B,, Bs, Bg, B10, B12» 12111211111111311111 3 32 880000 0.231 
Big, B23, Bro 1111111111 
8 B,, B;, Bg, B42, Baz, Bao 12111111111211221111 2 29 880000 0.242 
B33, B29 1111111111 
8 B4, Bs, Bg, Bio, Biz, Bag 12111111111211221111 2 29 880000 0.286 
B34, B29 1111111111 


4.3. Test system 3: New England 39-bus system 

New England 39-bus system is taken as test system 3. The information on the test system is shown 
in Table 7. The proposed algorithm is used to determine the minimum number and optimal locations of 
PMUs for the new England 39-bus system. In Table 8, the set with maximum redundancy value 51 and 
minimum is 49 ignoring ZIBs. After considering the setup cost of placement of PMU, the set with minimum 
redundancy value is considered for the best location for PMUs placement, ignoring ZIBs. Again, in Table 9, 
the installation cost of PMUs is considered. Therefore, the set with minimum redundancy value is considered 
for the OPP considering ZIBs. Bold letters show the optimal locations for the new England 39-bus system. 
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Table 7. System data of new England 39-bus system 
# of lines #of ZIBs _#of RBs Loc. of ZIBs Loc. of RBs 
46 12 9 B,, Bz, Bs, Bg, Bo, Bu, B30, B31, B32, B33, B34» 
B13, B14, B17, B19, B22 B35, B36, B37, Bag 


Table 8. Simulation results for the New England 39- bus system ignoring ZIBs 


# of PMUs Loc. of PMUs # of times each bus observed MR RA CBA($) Comp. time in s. 
12 By, By, Bg, Bo, B13, Bigs, 123121111111121 3 49 11,20,000 0.268 
Big B20, B22, Bos, Bag, 121212111231221 
Boo 111111111 
13 By, Bs, Bo, Bio, BaBa 121112111231211 3 51  11,80,000 0.249 
B17, B19, B20, B22, B23, 211221221121111 
B35, B29 111111111 
15 By, By, Bg, Bo, B13, Big, 113121111211121 3 51: 13,00,000 0.201 
Big B20, B22, Bos, Bre, 121312121211121 
B29, B32, B33, B36 111111111 


Table 9. Simulation results for the New England 39- bus system considering ZIBs 


# of PMUs Loc. of PMUs # of times each bus observed MR RA CBA($) Comp. time in s. 
8 B}, B4, B16, Bag, B20, Bos, 123211111111111 3 33 880000 0.231 
B36, B29 122211111231221 
111111111 
8 B}, Bg, Bio» By¢, B20 Bog, 121111111111111 2 33 880000 0.202 
B35, B29 111211112121111 
111111111 
10 By, B16 Big, B20, Bas, 112111111111111 3 35 10,00,000 0.228 
B26, B29, B32, B33, B36 121311111231221 
111111111 


4.4. Test system 4: IEEE 57-bus system 

IEEE 57-bus system is considered as test system 4. The data of the test system is given in Table 10. 
Using the proposed algorithm, the minimal number and optimal locations of PMUs for the IEEE 57-bus 
system are solved by ignoring and considering ZIBs. As shown in Table 11, the set with maximum and 
minimum redundancy values is 67 and 58. After considering the installation cost of PMUs, the set with 
minimum redundancy value is found suitable for the placement of PMUs ignoring ZIBs. The most promising 
locations for placement of PMU in IEEE 57-bus system ignoring ZIBs are shown in bold. In Table 12, the 
best PMU locations are {B,, Bo, B20, B29, B30, B32, B35, B44, Bso, B53, Bse} after considering the installation 
cost of PMUs. 


Table 10. System data of IEEE 57- bus test system 
#of lines #of ZIBs _#of RBs Loc. of ZIBs Loc. of RBs 
80 15 = By, B7, B11, B21, B22, B24, B26, B34, B36, B37, = 
B39, Bao, Bags Bag 


‘-- means not reported 


Table 11. Simulation results for the IEEE 57-bus system ignoring ZIBs 


# of PMUs Loc. of PMUs # of times each bus observed MR RA CBA($) Comp. time in s. 
15 B,, By, Bo, Boo, B24, 111111111111111 2 58 13,00,000 0.312 
By, B30, B32, Bas, Bao, 111111111211111 
B44, Bag, B50, B53, B56 211211111111111 
111111211112 
17 By, B3, B12, Byg, B21, B23, 122211111221112 2 61 14,20,000 0.352 
B34, Boo, B32, Bar, 111111222211111 
B33, Bas, Baz, Bsr, 211111111121211 
Bsa, Bs7 111111111121 
18 B,, Bs, B10, B14, B22, Bas, 111111111121213 3 67 14,80,000 0.361 
B37, B29, B30, Baz, Bas, 111111111211212 
B39, B53, B54 211212131121211 


111111222111 
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Table 12. Simulation result for the IEEE 57-bus system considering ZIBs 


# of PMUs Loc. Of PMUs # of times each bus observed MR RA CBA($) Comp. time in s. 
11 Bo, B20, B29, Bap, B32; 111111111111111 2 43 10,60,000 0.398 
B35, B44, Bso, Bs 3, B56 111111111111112 


112111111111111 
11111211111 


14 By, By, B12, Big, B23, Boo, 122211111221111 2 51 12,40,000 0.372 
B30, B32, B41, Baz, Bay, 111111111111111 
Bs, Bsa, B57 211111111121211 
111111111121 
15 B, Bs,Byo, By4, Bos, B27, 111111111121212 2 56 13,00,000 0.367 
B29, B30, B32, B35, B41, 111111111211212 
Baz, Bao B53, Bsa 211211111121211 


111111222111 


4.5. Test system 5: IEEE 118-bus system 

The data for the IEEE 118-bus test system is given in Table 13. The promising locations for the 
placement of PMU in the IEFE 118-bus test system ignoring ZIBs are shown in bold letters in Table 14. In 
this test system, the maximum redundancy is 4. Here, one bus is observed four times, three buses are 
observed three times, sixteen buses are observed two times and the remaining buses are observed once. 
Table 15 is showing the optimal locations of PMUs considering ZIBs. In this case, only twenty-eight PMU is 
required in order to monitor the IEEE 118- bus test systems. 


Table 13. System data of IEEE 118-bus system 
# of lines # of ZIBs #of RBs Loc. of ZIBs Loc. of RBs 
186 10 7 Bs, Bo, B30, B37, B38, Bos, Boa, Bog B71 Bio, B73, Baz, B111, B112, Bie Bury, 


Table 14. Simulation results for the IEEE 118-bus ignoring ZIBs 


# of PMUs Loc. of PMUs # of times each bus observed MR RA CBA($) Comp. time in s. 
32 By, Bs, Bo, B11, Byz, By7, Boy, 9 122221121121111211 4 148 23,20,000 0.698 
B24 B25, B28, B34, B37, 111111112111111211 
Bao, Bas, Bao 312212112111212112 


Bes, Bic, Be2, Bez, Beg Bz, 111121111112114211 

Bs, B37, Bso, Bos, Beg Boo, 112231122111211121 

111112111211211111 
1111111111 


Boa, Brow Bios, Biro, B114 


Table 15. Simulation results for the IEEE 118-bus system considering ZIBs 


# of PMUs Loc. of PMUs # of times each bus observed MR RA CBA($) Comp. time in s. 
28 By, Bg, By, By, Baz, B21:B25 121121111111112211 3 129 20,80,000 0.699 
Bzg, B33, B34, Bao, Bas, 111111113112111111 
B39, B52, B56» 311112112111212112 


Boo, B72, Bys, B37, Bgo, Bes, 111111111112113111 

Bse Boo, Boas Baoa, Baos,  112131121111211121 

111112111211211111 
1111111111 


Biro, B114 


The proposed algorithm is carried out successfully on different standard test systems, and the 
obtained results are compared with the various existing methods proposed in referred journals in Table 16 
(see Appendix). Based on the comparative analysis, it can be concluded that the proposed algorithm for OPP 
in the interconnected power network is more robust and efficient in terms of determining the minimum 
number and optimal locations of PMUs. 


5. CONCLUSION 

This paper presents a method for calculating the least number of PMUs to be installed and the 
appropriate locations for those PMUs to be set up in order to provide complete observability of the 
interconnected power network states under normal operating conditions. The A-star strategy, which is based 
on an admissible searching algorithm, has been developed for selecting the optimum site for the deployment 
of PMUs in order to provide comprehensive observability of the power network states. In the event that more 
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than one solution is available, a redundancy analysis is proposed, and the installation cost is taken into 
account when determining the most promising locations for the allocation of PMUs. A successful application 
of the suggested approach has been demonstrated on conventional IEEE 14-bus, IEEE 30-bus, New England 
39-bus, IEEE 57-bus, and IEEE 118-bus test systems, and the results have been compared to those achieved 
using other previously published approaches. The obtained findings demonstrate the efficacy and robustness 
of the suggested OPP technique for interconnected power network observability. 


APPENDIX 


Table 16. Comparative analysis of obtained results with standard published literature 


Test System 14-bus 30-bus 39-bus 57-bus 118-bus 
Chakrabarti and Kyriakides [4] 4 10 13 -- -- 
Roy et al. [5] 4 10 13 17 32 
Marin et al. [6] 3 T -- 12 29 
Venkatesh and Jain [8] 4 10 -- -- -- 
Tai et al. [9] 4 -- 13 -- 32 
Sodhi et al. [11] 4 -- 15 -- -- 
Mousavian and Feizollahi [12] 4 10 13 17 32 
Ramachandran and Bellarmine [13] 3 7 -- 12 28 
Gou [14] 4 -- -- -- -- 
Zhao et al. [15] -- 7 -- 11 29 
Liao et al. [16] 3 -- 8 11 28 
Sodhi et al. [17] 4 -- 15 -- -- 
Xia et al. [18] 4 10 -- 17 32 
Xu and Abur [22] 4 10 -- 17 32 
Xu and Abur [23] 3 -- -- 12 29 
Ahamadi et al. [26] 3 7 -- 13 29 
Proposed method 3 8 12 15 28 


‘--’ means not reported 
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